Comparing individual-based approaches to modelling the self-organization of multicellular tissues.

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage

BioSimulators: a central registry of simulation engines and services for recommending specific tools

Computational models have great potential to accelerate bioscience, bioengineering, and medicine. However, it remains challenging to reproduce and reuse simulations, in part, because the numerous formats and methods for simulating various subsystems and scales remain siloed by different software tools. For example, each tool must be executed through a distinct interface. To help investigators find and use simulation tools, we developed BioSimulators (https://biosimulators.org), a central registry of the capabilities of simulation tools and consistent Python, command-line and containerized interfaces to each version of each tool. The foundation of BioSimulators is standards, such as CellML, SBML, SED-ML and the COMBINE archive format, and validation tools for simulation projects and simulation tools that ensure these standards are used consistently. To help modelers find tools for particular projects, we have also used the registry to develop recommendation services. We anticipate that BioSimulators will help modelers exchange, reproduce, and combine simulations

Early Embryonic Vascular Patterning by Matrix-Mediated Paracrine Signalling: A Mathematical Model Study

During embryonic vasculogenesis, endothelial precursor cells of mesodermal origin known as angioblasts assemble into a characteristic network pattern. Although a considerable amount of markers and signals involved in this process have been identified, the mechanisms underlying the coalescence of angioblasts into this reticular pattern remain unclear. Various recent studies hypothesize that autocrine regulation of the chemoattractant vascular endothelial growth factor (VEGF) is responsible for the formation of vascular networks in vitro. However, the autocrine regulation hypothesis does not fit well with reported data on in vivo early vascular development. In this study, we propose a mathematical model based on the alternative assumption that endodermal VEGF signalling activity, having a paracrine effect on adjacent angioblasts, is mediated by its binding to the extracellular matrix (ECM). Detailed morphometric analysis of simulated networks and images obtained from in vivo quail embryos reveals the model mimics the vascular patterns with high accuracy. These results show that paracrine signalling can result in the formation of fine-grained cellular networks when mediated by angioblast-produced ECM. This lends additional support to the theory that patterning during early vascular development in the vertebrate embryo is regulated by paracrine signalling

A Novel Method to Verify Multilevel Computational Models of Biological Systems Using Multiscale Spatio-Temporal Meta Model Checking

Insights gained from multilevel computational models of biological systems can be translated into real-life applications only if the model correctness has been verified first. One of the most frequently employed in silico techniques for computational model verification is model checking. Traditional model checking approaches only consider the evolution of numeric values, such as concentrations, over time and are appropriate for computational models of small scale systems (e.g. intracellular networks). However for gaining a systems level understanding of how biological organisms function it is essential to consider more complex large scale biological systems (e.g. organs). Verifying computational models of such systems requires capturing both how numeric values and properties of (emergent) spatial structures (e.g. area of multicellular population) change over time and across multiple levels of organization, which are not considered by existing model checking approaches. To address this limitation we have developed a novel approximate probabilistic multiscale spatio-temporal meta model checking methodology for verifying multilevel computational models relative to specifications describing the desired/expected system behaviour. The methodology is generic and supports computational models encoded using various high-level modelling formalisms because it is defined relative to time series data and not the models used to generate it. In addition, the methodology can be automatically adapted to case study specific types of spatial structures and properties using the spatio-temporal meta model checking concept. To automate the computational model verification process we have implemented the model checking approach in the software tool Mule (http://mule.modelchecking.org). Its applicability is illustrated against four systems biology computational models previously published in the literature encoding the rat cardiovascular system dynamics, the uterine contractions of labour, the Xenopus laevis cell cycle and the acute inflammation of the gut and lung. Our methodology and software will enable computational biologists to efficiently develop reliable multilevel computational models of biological systems

A new mechanism for collective migration in Myxococcus xanthus

Myxobacteria exhibit a complex life cycle characterized by a sequence of cell patterns that culminate in the formation of three-dimensional fruiting bodies. This paper provides indications that the specific cell shape of myxobacteria might play an important role in the different morphogenetic processes during the life cycle. We introduce a new mechanism for collective migration that can explain the formation of aligned cell clusters in myxobacteria. This mechanism does not depend on cell cooperation, and in particular it does not depend on diffusive signals guiding cell motion. A Cellular Potts Model (CPM) that captures the rod cell shape, cell stiffness and active motion of myxobacteria is presented. By means of numerical simulations of model cell populations where cells interact via volume exclusion, we provide evidence of a purely mechanical mechanism for collective migration, which is controlled by the cells' length-to-width aspect ratio

Influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics within human airway epithelium.

Human airway epithelium (HAE) represents the primary site of viral infection for SARS-CoV-2. Comprising different cell populations, a lot of research has been aimed at deciphering the major cell types and infection dynamics that determine disease progression and severity. However, the cell type-specific replication kinetics, as well as the contribution of cellular composition of the respiratory epithelium to infection and pathology are still not fully understood. Although experimental advances, including Air-liquid interface (ALI) cultures of reconstituted pseudostratified HAE, as well as lung organoid systems, allow the observation of infection dynamics under physiological conditions in unprecedented level of detail, disentangling and quantifying the contribution of individual processes and cells to these dynamics remains challenging. Here, we present how a combination of experimental data and mathematical modelling can be used to infer and address the influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics. Using a stepwise approach that integrates various experimental data on HAE culture systems with regard to tissue differentiation and infection dynamics, we develop an individual cell-based model that enables investigation of infection and regeneration dynamics within pseudostratified HAE. In addition, we present a novel method to quantify tissue integrity based on image data related to the standard measures of transepithelial electrical resistance measurements. Our analysis provides a first aim of quantitatively assessing cell type specific infection kinetics and shows how tissue composition and changes in regeneration capacity, as e.g. in smokers, can influence disease progression and pathology. Furthermore, we identified key measurements that still need to be assessed in order to improve inference of cell type specific infection kinetics and disease progression. Our approach provides a method that, in combination with additional experimental data, can be used to disentangle the complex dynamics of viral infection and immunity within human airway epithelial culture systems

Integrating intracellular dynamics using CompuCell3D and bionetsolver:Applications to multiscale modelling of cancer cell growth and invasion

In this paper we present a multiscale, individual-based simulation environment that integrates CompuCell3D for lattice-based modelling on the cellular level and Bionetsolver for intracellular modelling. CompuCell3D or CC3D provides an implementation of the lattice-based Cellular Potts Model or CPM (also known as the Glazier-Graner-Hogeweg or GGH model) and a Monte Carlo method based on the metropolis algorithm for system evolution. The integration of CC3D for cellular systems with Bionetsolver for subcellular systems enables us to develop a multiscale mathematical model and to study the evolution of cell behaviour due to the dynamics inside of the cells, capturing aspects of cell behaviour and interaction that is not possible using continuum approaches. We then apply this multiscale modelling technique to a model of cancer growth and invasion, based on a previously published model of Ramis-Conde et al. (2008) where individual cell behaviour is driven by a molecular network describing the dynamics of E-cadherin and -catenin. In this model, which we refer to as the centre-based model, an alternative individual-based modelling technique was used, namely, a lattice-free approach. In many respects, the GGH or CPM methodology and the approach of the centre-based model have the same overall goal, that is to mimic behaviours and interactions of biological cells. Although the mathematical foundations and computational implementations of the two approaches are very different, the results of the presented simulations are compatible with each other, suggesting that by using individual-based approaches we can formulate a natural way of describing complex multi-cell, multiscale models. The ability to easily reproduce results of one modelling approach using an alternative approach is also essential from a model cross-validation standpoint and also helps to identify any modelling artefacts specific to a given computational approach